o 



Time-dependent Ginzburg-Landau Equation in the 
Nambu-Jona-Lasinio Model 

Wei-jie Fu'^, Da Huang^, Fa-bing Wang*^ 

""Kavli Institute for Theoretical Physics China (KITPC), Key Laboratory of Frontiers in 
fvq ' Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Science, 

►^ . Beijing 100190, China 

o 

:z; 

o 

<^ ■ Abstract 

r-| ! We apply the closed time-path Green function formalism in the Nambu- 

Ph| Jona-Lasinio model. First of all, we use this formalism to obtain the well- 

jil^' known gap equation for the quark condensate in a stationary homogeneous 

system. We have also used this formalism to obtain the Ginzburg-Landau 

(GL) equation and the time- dependent Ginzburg-Landau (TDGL) equation 

^ \ for the chiral order parameter in an inhomogeneous system. In our derived 

^ \ GL and TDGL equations, there is no other parameters except for those in 

the original NJL model. 
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1. Introduction 

QCD thermodynamics, for example the equation of state of quark gluon 
plasma (QGP), phase transition of chiral symmetry restoration, deconfine- 
ment phase transition and so on, has been a subject of intensive investiga- 
tions in recent years. On the one hand, the deconfined QGP are expected 
to be formed in ultrarelativistic heavy-ion collisions [l|, |2|, y, y, |5|, la, LZl, |8[ 
(for example the experiments at the Relativistic Heavy Ion Collider (RHIC) 
and at the Large Hadron Collider (LHC)) and in the interior of neutron 
stars [qI, llO, lul, ll2|; On the other hand, studying the thermodynamical be- 



haviors of the QGP, especially the deconfinement and chiral phase transitions, 
is an elementary problem in strong interaction physics. 
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Lattice QCD calculations [l3|, llj, ll5|, [l6| as well as models studies [l7|, 



18l . Il9l . I20I . I21I . I22I . I23I . I24J . |25| indicate that there is a critical point in the 
QCD phase diagram in the plane of temperature and baryon chemical po- 
tential 261], which separates the first order phase transition at high baryon 
chemical potential from the continuous crossover at high temperature. Fur- 
thermore, experiments with the goal to search for the QCD critical point are 
planed and underway at RHIC and at the Super Proton Synchrotron (SPS) 



27|, |28|, l29|, |30 



Therefore, the critical dynamics of the chiral phase transi- 
tion has attracted lots of attentions in recent years. Based on equilibrium 
thermodynamics 3l| or quasi-stationary framework [32i|, signatures of the 
QCD critical point have been studied by M. A. Stephanov et al.. A. Singh et 
al. have also studied the kinetics of the chiral phase transition subsequent 
to a quench from a disordered phase to the ordered phase 331 . 

In this work, we will adopt the closed time-path Green function (CTPGF) 
formalism to derive the Ginzburg-Landau equation and the time-dependent 
Ginzburg-Landau equation of the order parameter for the chiral phase transi- 
tion in an inhomogeneous quark matter system in the Nambu-Jona-Lasinio 
(NJL) model. The CTPGF formalism, developed by Schwinger [341 and 



The CTPGF formalism, developed by Schwinger 
Keldysh [35[, has been used to solve lots of interesting problems in statistical 
physics and condensed matter theory ^], and it has also been used in the 
NJL model to derive the transport equations [STJ. It is generally believed that 
this technique is quite effective in investigating the nonequilibrium statistical 



theory [36|, |38(]. It has also been used to treat a system of self-interacting 
bosons described by Xcj)^ scalar fields [391] . 

The paper is organized as follows. In Sec. [2] we simply review the CTPGF 
formalism, mainly on the generating functional of the CTPGFs. In Sec. [3] we 
apply the CTPGF formalism into the NJL model. First of all, we use this 
formalism to derive the well-known gap equation of the quark condensate for 
the homogeneous system, then we also obtain the Ginzburg-Landau equation 
and time-dependent Ginzburg-Landau equation of the order parameter for an 
inhomogeneous system. In Sec. H] we present our summary and conclusions. 



2. Simple Review on the CTPGF 

2.1. Two-point CTPGFs 

In this section we give a short review about the CTPGF formalism and in- 
troduce some notations which we will use in our following discussions. Read- 
ers who are not familiar with the CTPGF formalism are strongly suggested 



to reference the excellent review article by K. C. Cliou et al. 36|. Here, we 



take the real boson field ^{x) for example. The two-point CTPGF is defined 

as 

G{x,y) = -zTi {Tp{ip{x)ip{y))p} 

^ -z{T,{v{xMy))), (1) 

where p is the density matrix and Tp is the time ordering operator along the 
closed time-path p which goes from — oo to -|-oo (also denoted as the positive 
time branch) and then returns back from -|-oo to — oo (negative time branch). 
We should note that any event at the negative time branch is later than any 
that at the positive time branch. The Green function G{x,y) in Eq.(IT]) on 
the closed time-path p can be expressed as Green function whose coordinates 
X and y are constrained to be on either positive or negative time branches, 
and then G{x,y) becomes a 2 x 2 matrix, i.e., 

Gix,y)^(^-- ^^-1-f^" ±) (2) 



^ G-+ G— J ~ \ G- G^ 

with 

GF{x,y) ^ -t{T{^ixMy))), (3) 

G+{x,y) = -i{if{y)if{x)), (4) 

G-{x,y) = -i{^{x)^{y)), (5) 

Gp{x,y) ^ -t{f{v{xMy))), (6) 

where T is the usual time-ordering operator and T is anti-time-ordering op- 
erator. Therefore, Gp is the conventional Feynman causal propagator and 
Gp is an anti-causal propagator, whose expression can be given explicitly as 

Gpix, y) = -i9{yo, Xo){(pix)ip{y)) - ie{xo, yo){(piy)(pix)). (7) 

Furthermore, It can be easily proved that we have the following identity, i.e., 

Gp{x, y) + Gp{x, y) = G+{x, y) + G_(a:, y). (8) 

So far we have presented CTPGFs in two representations: one is the closed 
time-path representation, i.e., Eq.([T]) and the other is the single time repre- 
sentation given by Eq.(|2]). In fact, there is another representation which is 



directly related with measurable quantities, which is often called as physical 
representation and defined as 

Gr{x,y) = -ie{xo,yo){[(p{x),ip{y)]), (9) 

Gaix,y) = ie{yo,xo){[(pix),ip{y)]), (10) 

G,ix,y) = -i{Mx),^iy)}). (11) 

Or, Ga, and Gc are retarded, advanced, and correlation Green functions, re- 
spectively. The relations between the CTPGFs in the physical representation 
and those in the single time representation are given by 

Gtj. = (jr p — G-\- = Gr_ — Gp, [.^■^) 

Ga = Gp — G- = G^ — Gp, (13) 

G, = Gp + Gp = G+ + G^, (14) 

and the inverse relations are 

Introducing two-component vectors 
ones can express Eq.( !T5|) as 

Gal3 = -GriaVli + -Gaf]a^l3 + -Gc^a^p, (17) 

where Gap with Greek subscripts a, /3 = ± are components of the CTPGF 
matrix in single time representation in Eq.([2]). 

2.2. Generating junctionals of CTPGFs 

We begin with the Lagrangian density of the real boson field as follows 

C = -d^ip{x)d^(p{x) - -w?(p^{x) + Cint{<^{x)) 

= Co{^{x)) + C,nt{^{x)), (18) 



where Cq is the free field term and £j„j is the interaction term. The generating 
functional for the CTPGFs is defined as 



Z[J{x)\ = Tr< Tp exp ii I d'^xJ{x)(p{ 
where the time integration is 



(19) 



/+00 n—oo 

dt+ + / dt. 
,_ 'Oo J +00 

+00 



dt+- 



dt. 



(20) 



The n-point CTPGF is defined as 



G,{l---n) = (-z)"-^Tr 



^(-1)' 



Tp{ip{l)---ip{n))p 
6'^Z[J{x)] 



6J{l)---6J{n) 



j=o 



(21) 



In the interaction picture the generating functional Z[J(x)] in Eq.( IT9ll can 
be expressed as 

Z[J{x)] = TrJTp exp \^i / d'^x{Cint{^Pi{x)) + J{x)ipj{x))j pj 



exp 



i I d^xCr 



int 



6J{x) 



X Tr<^ T, 



exp [id xJ{x)(pi{x) 
p 



(22) 



Using the Generalized Wick theorem J40| , one can show [4l|, |42|, |36 



T 

■J- n 



with 



expli / d'^xJ{x)(pj{x)] = Zo[J{x)] : exp i / d'^xJ{x)(pi{x) :, (23) 



(24) 



Zo[J{x)] = / [d(p{x)]exp id x[CQ{ip{x)) + J{x)(p{x) 
J p J p 



Here : : denotes the normal product. Substituting Eg. ( 12^ into Eg. ( 12^ . one 
obtains 



Z[J(x)] = exp 



I if^^^nt[-ljj^^ 



Z,[J{x)]N[J{x)], (25) 



here 



N[J{x)\ = Tr< : exp i I d'^xJ{x)(pi{x) : f)> 



^ e^^{iW^'[J{x)]) 



(26) 



with 



°° -n-l 

6 



n=l ^ 



dl--- / dnJ{l) ■ ■ ■ J(n)Tr[: </?/(!) ■ ■ ■ ifii{n): p]. 



(27) 
where A^[J(x)] is the correlation functional for the initial state. 

After some calculations, the generating functional Z[J(x)] in Eg. (12 5 p can 



be expressed as [42 



Z[Jix)] 



[d(^ix)] exp i / d^x(Co{(f{x)) + Cint{^{x)) + J{x)ip{x] 



X exp 



iW, 



N 



d vGq {x,y)ip{y) 



(28) 



where Gq ^ is the inverse of the propagator of the real boson field, given by 
Go\x,y) = i-d^d^ - m^%{x - y). (29) 

In the single time presentation 5p{x — y) can be expressed as 

5p{x -y) = 5{x - y)a^. (30) 

The generating functional for the connected CTPGFs is given by 

W[J{x)] = -i\nZ[J{x)]. (31) 

Correspondingly, one can obtain 



Gl{l---n) = (-z)"-^Tr 



Tp{ip{l)---ip{n))p 
1 6"W[J{x)] 



6J{l)---6J{n) 



j=o 



(32) 
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For n = 1 we have 

(.(.)) ^ ^™. (33) 

dJ[x) 

The vertex generating functional can be obtained by performing the Legendre 
transformation as follows 



r[((^(x))] = W[J{x)] - / d''xJ{x){<^{x)). (34) 

Jp 

It can be easily found that 

S{<^{x)) 

Next, we express the external source term in the generating functional in 
Eg. (128!) in the single time representation as 

d'^xJ{x)ip{x) 

dtSx{J^{x)Lp^{x) — J^(x)ip-{x)) 

oo 

(i^xJ^(x)(T3(^(x), (36) 

with 



p 

+ 00 



#wH^:{:|), ^w-(':g). (3^) 



Introducing 



J^{x) = J+{x) - J-{x), Jc{x) = -{J+{x) + J_{x)), 
^a{x) = ^+{x)-ip-{x), ipc{x) = -{ip+{x) + ip^{x)), (38) 
one can express Ig as 

Is = d'^x{Ji,{x)(pc{x) + Jc{x)(p^{x)). (39) 

Then we can obtain the generating functional for the CTPGFs in the physical 
presentation as follows 

Z[J^{x), Je(x)] = Z[J4x), J4x)]. (40) 



In the same way, we have 

W[J^{x),Ux)] = -i\YiZ[J^{x),J,{x)]. (41) 

and 

V[{ip^{x)),{ip,{x))] = W[J^{x),Ux)] 

d^xi.Mx)Mx)) + Jc{x){ip^{x))), (42) 



with 



i^cix)) 



5W[.Mx),.Ux)] 



5J^{x) 
Therefore, one can also obtain 

5T[{ip^{x)),{ipc{x))] 



(^,(,)) = mJ.(^),M^)]^ (43) 



5{'^^{x)) 



Jc{x), 



6Jc{x) 



5T[{ip^{x)),{ipcix))] 



-Mx). (44) 



3. Application of the CTPGFs in the NJL model 



In this section we use the CTPGFs to study the spacially and temporally 
inhomogeneous quark matter in the NJL model. The Lagrangian density for 



the two-flavor NJL model is given by [43|, |44, |45|, |46|, |47|, |48 

^NJL 



i){i-i^d^ - mQ)i) + G {^jj^f + {ipi-i^Tipf 



(45) 



where ip = {ipu,4'd)'^ is the quark field, and mo = diag(m„,m(i) is the cur- 
rent quark mass matrix. Throughout this work, we take rriu = rrid = mo, 
assuming that the isospin symmetry is reserved on the Lagrangian level. The 
four-fermion interaction with an effective coupling strength G for the scalar 
and pseudoscalar channels has SUv(2) x SUa(2) x Uv(1) symmetry, which 
is broken to SUv(2) x Uv(l) when mo ^ 0. Here T'^{a = 1,2,3) in the 
Lagrangian density in Eq. fj45|) are Pauli matrices in flavor space. 

Following the CTPGF formalism in Sec. [21 we can obtain the generating 
functional for the NJL model as follows 

Z[h,h'';f],r]] = / [dV'][o?V'] exp jz / c?^x(^£o(V^, V^) + Ant(V^,^) + # + ^^7 



+ /i^V^ + /,>75r> + W^j^[-V^1So\-^oV]) 



(46) 



with 






ij{ijf,d^ - rho)ip, 



G 



(ijipf + {ipi'j^Tipy 



(47) 



where f], t] are the external sources for field ip and ip, respectively; h, /i" for 
composite operators ijjip and tpi'-^^T'^ip . The term W^ in Eq. fH6|) is due to 
initial correlations as shown in the last section, and here S^^ is the inverse 
of the propagator for the free fermion field. Up to a constant, we have 



exp \i / d x\G 



(ipipf + {i)i'j5Tipf + h^jip + h''ipi'j5T''ip 



[da] [dn"'] exp <i d'^x 
p •>p 



G{a 



2 I ^a.2\ 



{2Gip'4j + h)a 



+ (2GV^75r> + /i'')7r" 



4G 



{h^ + h 



a2' 



(48) 



Substituting Eq.l HHj) into Eq.( H6|) . we obtain 



-y^.'.(ft^ + ft-)i 



[daWd-K"] I [diPWdiP] 

V 



Z[h,h"';ri,ri] = exp 

X exp h I d^x (t/;^- V - G{a^ + vr"^) + ha + h^ir" 

+ r]i; + i;r^ + W^''[-ijX\-X'i^])}, (49) 



with 



S~\x, y) = [i^^d^ - (mo - 2Ga) + 2Gi^^,ti''t''%{x - y). (50) 

In the realistic calculations for the equilibrium (or local equilibrium) system, 
the contribution of the density matrix, i.e., the term W^ in Eq. (l49|) can be 
neglected temporarily for simplicity, and its contribution can be recovered 
through the fluctuation-dissipation theorem satisfied by the CTPGFs. More 
detailed discussions can be found in Refs. [36|,|49|. Neglecting W^ in Eq. (jl9|) 
temporarily, we can integrate over the fermion field and obtain 

Z[h, /i^ f],r]] = exp [ - ^ /" d^xih^ + /i"')] f[da] [dn"] 

/o(a,7r'^) + Jdet(a,7r») 



X exp < t 

4 



+ 



d''x{ha + /iV^) - / d''xd^yf]{x)S{x,y)r]{y) 
p Jp 



(51) 
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with 

/o(a,7r") = - fd^xGia^ + n"^), 
Jp 

/dct(a,7r") = -iTTlnS-\ (52) 

where Tr means trace operation in both the inner space (Dirac, flavor, and 
color) and the coordinate space. Then we can introduce the effective action 
for the a and vr flelds as follows 

P^'ia, tt'^) = /o((T, tt") + /det(a, tt"^). (53) 

In the mean fleld approximation, the vertex generating functional can be 
approximated by the effective action '50|], i.e., 

r(cr,7r")~r^(a,7r"), (54) 

where a and n"- mean (a) and (tt"), without confusions. In the physical 
representation we have 

r(o-A,<;crc,0 ~ r'^(crA,7r";(Tc,7r°) 

= /o(aA,<;ae,7r;f) + /det(aA,<;ae,0, (55) 



with 



and 



Jo((TA,<;a„0 = -2G /"rf^a:(aAa, + «), (56) 



+ 2Gz75r'^(ea< + lvaK)v^] SapS\x - y). (57) 

Following the notations in the last section, we have 

a^{x) = cr+(x) - cr_{x), a^x) = -(cr+(x) + a_(x)), 

<(x) = 7r^(x)-7r^(x), <(a;) = ^«(a;) + 7r^(x)). (58) 

In the following, we neglect the influence of the vr" fleld and concentrate on 
the order parameter of chiral phase transition, i.e., the a fleld. First of all, 

10 



let us consider the simplest stationary homogeneous quark system, i.e., the 
a field is independent of the time and space coordinates. From Eq. (144 p and 
Eq.f lM|) we have 



Sa^ 



-K = 0, 



cr^=0,ac=cr 



For the left hand side of Eq. (159!) , we have 









6a^ 






SIdct{(^A,0-c) 








6a^ 


Eq.iZ), 


we have 






6[S- 


\y, 


z)U _ ns 



IGo,, 



-iYy[S 



.6S 



-1 



6aA 



SaJx) 



G5^pS\y-x)S\y-z). 



(59) 



(60) 



(61) 



(62) 



Then 



SaA(x) 



—i / d'^yd'^z tr S{y,z 
Jp 

—i / d'^yd'^z ti 
—i / d'^yd'^z ti 



5S-'{y,z) 



alS{y, z)a^ 
S{y,z 



S(Ta{x) 
p6S-\y,z) 



6ai:,{x 
6S-\y,z) 



6(7. (x) 



(63) 



where tr means trace operation not including in the coordinate space; a^ is 
the Pauli matrix in the single time representation. Substituting Eq.( !T7|) and 
Eq.( l62|) into above equation, we obtain 



6Idet{(^A,(^c) 

6(7. (x) 



-iGtiSJx, x). 



(64) 



For the inverse of the quark propagator in Eq. (l50|) without vr" field, we can 
easily obtain 



Srip) - Sa{p) = S^{p) - S+{p) 



-ivr^ — [6{po - Ep) - 6{po + Ep)], 



E., 



(65) 



p 
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in the inomentuin space. Here we have 



M = mo - 2Ga, Ep = ^/fTM^. 



(66) 



Employing the fluctuation-dissipation theorem [36[, we obtain 



Sc{p) 



tanh [-(po - IJ^)][Srip) - Sa 



-m [(5(po 



Er, 



SiPo 



Ep)] tanh [-{po 



/i)], (67) 



where P = 1/T is the inverse of the temperature, fi the quark chemical 
potential. Substituting Eq.f l67|) into Eq. flMj) . we obtain 

5Idet{crA,(^c) ,^ f d'^P 



Scr^{x) 



-%G 



(2vr) 



-\GNjN^ 



ftr5c(p), 

S-p M 
WfYp 



1 



^P(E^-ii) _^ ]_ 



where N f is the number of quark flavor, N^. the color number. Combining 
Eqs.([59D(l60D with Eq.(l6HD, we obtain 



a 



-2NfN, 



d^p M 



1- 



1 



oJ3{Ep-ij) _|_ I 



(69) 



which is the gap equation for the order parameter in the NJL model. 

In the following we will consider the weak non-stationary inhomogeneous 
quark matter, where the "weak" means that we can expand the left hand side 
of Eq. flMJl) as power series of the gradient of the a{t, x) fleld. The efficiency of 
this expansion has be verified in Ref. 49|]. Furthermore, we only concentrate 
on the quark system which is near the chiral phase transition, therefore, the 
order parameter, i.e., the expectation value of the a{t, x) field is smaller than 
the critical temperature. Then, the left hand side of Eq. ( I59l) can be expressed 
as 



51 



O'AjO'c 



6a Ax) 



(^tl^CFc 



O"A=0,CTc=(t(x) 



+ 



5a Ax) 



dS 



= 0,fJc = 



(52/'=«(a„a,) 



5(7/ 



{x+ \)5ac{x- \] 



a[x 



2' 






(70) 
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where the dependence of the a field in the fermion propagator on the coor- 
dinate is neglected in the first term of the right hand side of Eg. (170]) . i.e., it 
is just the result for the homogeneous system with a replaced by a{x). The 
second term of the right hand side of Eg. (1701) can be expressed as 



+ 



1 



d'y 
d'y 



S'P'^icT^.CT,) 



6a^{x + ^)S(Jc{x - 


I) 


52/^ff(a„a,) 




6ai^{x + ^)S(7c{x - 


1) 


S'P^'ia^^a,) 








y'd,a{x) + -y>'y''d^d,a{x) 



6ai^{x + \)6(Jc{x - \] 



— I 



_d_ 



o^iy 



q=0 



d^a{x) 



92 



o^iy 



2' "' ydq^dq^ 
d 



q=0 



df,di,a{x) 



Idq, 



-Tr{x,q) 



q=0 



df.(r{x) - - 



d^ 



dqndqy 



Tr{x,q) 



q=0 



d^du(T{x) 



with 



Tr{x,q) = / d'^y 



s'r^'i^^^c) 



-JQV 



6ai,{x + ^)5ac{x - |) 



(71) 



(72) 



where /x, u =0, 1, 2, 3. We have 



+ 



5^/det(0"A,<7c) 



5ac,{x + l)6ac{x-l) 5a^{x + l)6a^{x - D ^^^(a; + f )(5(Tc(x - |) 
with 



5 Io{ai,,ac 



and 



(5(Ja(x + ^)6(7c{x - I 



2G6\y), 



(73) 
(74) 



6ai,{x + ^)Sac{x- |) 
After some calculations, we find 



zTr 



S 



S- 



6S-^ 



Sai,{x + ^) 5ac{x-^) 



(75) 



Tr{x,q) = -2G + t2G^ 



d^p 



(2vr; 



^ti[Scix,p)Sr{x,p + q) + Saix,p)Sc{x,p + q)]. 

(76) 
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The correlation Green function for quark field Sc is given in Eq.f l67p . and the 
retarded and the advanced Green functions are given by 



1 i^ + M 



^- M + ie E. 



V)i 



V 



Here, the dependence of the Green functions on the space and time coor- 
dinates is realized through the (y{x) field as shown in Eq. fl66|) . Performing 
calculations on Eq. flT^ . we obtain 

Tr{x,q) = -2G 

1 f 2M^ + Epqo-p-q 



Ep V (go + E^y - E^^g + ie{qo + Ep) 
, 2M'-Epqo-p-q x 

(q,-Epy-E'p_,^ + te{qo-Ep)J' ^ ^ 

When considering the variation with respect to the space coordinates, we can 
simplify Eq. fl78|) further by taking the limit qo ^ firstly, which gives 



+ AG'N^Nf / ^(tanh [^{Ep - ^^)] + tanh [^{Ep + ^^)] 
2M^-p-q/ 1 ^ 1 \ 

^p^p+q ^ ^p ~ ^p+q ^p "T -C/p+g ^ 

where V denotes main value, and 

Ep^q=[{p + q)' + M'Yl\ (80) 

when |g| — )■ 0, we can expand the left hand side of Eq. (179|l as power series 
of |g|. Furthermore, taking into consideration that the system is near the 
chiral phase transition, so the quark mass M is smaller than the critical 
temperature, we obtain 

Tr{x,q)q^=o ~ -20 
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+ AG'N.Nf— / dp^ 



vr^ 



A p2 



Ep L e/5(-^p-^) + 1 eP^'^v+t^) + l 



IgrC^iV.iV;- rfp- 



vr^ 



^« 



(81) 



where A is the cut-off of the momentum integration which is a parameter in 
the NJL model. From Eq. flHT]) we can easily obtain 



|-r.(x,g) 



0, 



92 



dq^dq- 



-Tr{x,q) 



9=0 



9=0 



(82) 



1 



-S^,^2G'N,Nf- I dp— 



■n" 



1 



1- 



,P{Ep^ii) j^ I 



1 



g/3(£;p+M) _^ 1 



(83) 



Substituting Eqs.(IH2l)(IH3l) into Eqs. ([70]) (JTID and combining Eq.(l59]), we ob- 
tain 



1 f^ M 
a{x) + N,Nf— / dp/ — 

71" Jo ^p 


-iT^cNf-l dp- 

2 TT^ Jq hp 


1- 


= 0, 





1 1 



^l3{Ep-tM) _^ I ^p{Ep+fi) _^ I 



VV(x) 



(84) 



which is the Ginzburg-Landau equation for the inhomogeneous quark matter 
in the NJL model. 

We can also consider the time variation of the quark matter near the criti- 
cal point by the same method that we applied to the space coordinates, which 
gives us the time-dependent Ginzburg-Landau equation. Without influenc- 
ing the final result, we can take the limit if — ;■ and Eq. (l78p is simplified to 
the following form. 



Tr{x,q)q-=Q = -2G 

+ AG'N^Nj J (1^ ( tanh [^{Ep - /.)] + tanh [^{Ep + /i)] 
2M2 + Epqo 



d^p 

2M2 
Ep V (go + EpY^El + 2e(go + Ep) 
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+ 



2M2 - EpQo 



[q,-E,f-El + tt{q,-E,] 



(85) 



In the same way, we ignore the mass term because the system is near the 
QCD critical point. By expanding Eq.(l85l) with respect to go, we obtain 



r^(x,g)g-=o 



-2G 



1 r^ . p2 



-K^ In E, 



1 







1- 



1 1 



^l3{Ep-tM) _^ I ^l3{Ep+n) -\-l_' 

(86) 



Differentiating Tr{x,q) with respect to go, we have 
d 



Idq' 



o'^rix.q) 



Jg=0 



dq^dq^ 



-Tr{x,q) 



q=0 



0, 



2G^NcNf 



(87) 



1 r"" . p^ 



El 



Tl" 



dp 



1- 



1 



Q(3{Ep-fi) _^ I 



g/3(Bp+M) _^ I 



(88) 



Substituting Eqs.(IH7])(IHHD into Eqs.(I70]) (IZID and combining Eq.(l8l]), we ob- 
tain the time- dependent Ginzburg-Landau equation: 



aix) + N,Nf — 



TX^ 



—rN.Nf — 



G 
~2 
G 
2" 
0. 



dp 







dpp' 



E„ 



+^N,N,- 



TT 



1 f\ p' 
-^p 



M 
Ep 

1- 
1- 



1 



1 



1 1 



1 1 



Ql3{Ep-ii) _^ I ^p{Ep+ii) _j_ I 



VV(s) 
dla{x) 



(89) 



Note that in the above time-dependent Ginzburg-Landau equation, the time 
derivative is second order, meaning that the a mode is a propagating one, 
which is different from the diffusive time-dependent Ginzburg-Landau equa- 
tion for the precursory mode in the superconductor phenomena 5l|, l52 . 
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4. Summary and Discussions 

CTPGF formalism is a powerful tool, because it can deal with not only 
equilibrium problems but also nonequilibrium ones. When much attention 
are paid on the critical dynamics of the QCD phase transitions, especially the 
chiral phase transition, it is quite natural for people to choose the CTPGF 
formalism. In this work, the CTPGF formalism is applied in the NJL model. 
First of all, we simply review the CTPGF formalism, mainly on the gener- 
ating functional of the CTPGFs. Then we use this formalism to obtain the 
well-known gap equation for the quark condensate in stationary homogeneous 
system. We have also used this formalism to obtain the Ginzburg-Landau 
equation for the chiral order parameter in stationary inhomogeneous system 
and the time- dependent Ginzburg-Landau equation describing the critical 
dynamics of the chiral phase transition. In our derived Ginzburg-Landau 
equation, there is no other parameters except for those in the original NJL 
model. 

In the present paper, we only derive the time-dependent Ginzburg-Landau 
equation. However, to fully appreciate the critical dynamics of the chi- 
ral phase transition, we also need to analyze the equation numerically and 
make some comparison with the experimental data. Furthermore, the present 
derivation is only the mean field result and we ignore the effects of the in- 
teraction of (t{x) mode with 7r(a;) meson. Related works are under progress 
and we will report them elsewhere. 
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